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Abstract 



I In this paper we present a study of magnetic order of quantum Heisenberg spins 

Pi ■ on pyrochlore type lattices: a lattice of corner sharing tetrahedra. We use the Cluster 

Variation Approximation with the tetrahedron as basic cluster. We find that in case 
of anti-ferromagnetic nearest-neighbour exchange there is no long-range spin order 
at any finite temperature for any spin s. For ferromagnetic exchange and s = 1/2 
^> . the spin system remains in a non-magnetic phase for all temperatures, while for 

0\\ s > 1/2 there is a finite order temperature. We derive the ground state phase 

diagram as a function of exchange anisotropy and tetragonal anisotropy and find 
OO ■ regions of re-entrant behaviour. Magnetic fields can induce one or more magnetic 

, transitions. Schottky-like anomalies appear in the specific heat, with on top signs 

I of re-entrant behaviour. The limits of this Cluster Variation Approximations are 

discussed. 
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1 INTRODUCTION. 



In this paper we report results for the critical temperature, susceptibility, magnetiza- 
tion and specific heat for Quantum Heisenberg spins (s = 1/2 to s = 7/2) on the highly 
frustrated pyrochlore lattice using the Cluster Variation Approximation (CVA). This py- 
rochlore lattice is a lattice with Corner Sharing Tetrahedra (CST) which has the highest 
known degree of frustration according to Lacorre's definition [|l|. 

In the last decade the properties of spins on this pyrochlore lattice (see figure |1|) have 
received much attention because it promises to be a realization of a 3D Quantum Spin 
Liquid (QSL) with possibly exotic ground states. Investigated Quantum ground states 
are for example based on singlet valence bonds, leading to a Resonating Valence Bond 
(RVB) system, introduced by Anderson |^ for a 2D triangular lattice, or a Valence 
Bond Solid/Crystal, or order based on a singlet tetrahedron (=RVT) or plaquette state, 
or the formation of a ground state with long range chiral order. 

One of the basic features of such frustrated spin lattices is their extensive residual or 
zero-point entropy at T = 0, or formulated in another way - a macroscopic degenerate 
ground state, with a degeneracy oc a^, where a > and is of the order of Avogadro's 
number. The consequences of this degeneracy are many: 1) a large decrease of the ordering 
temperature Tc, such that / = Qcw /Tc S> 1, where Qcw is the asymptotic Curie- Weiss 
temperature [|, in the region below the classical Mean-Field transition temperature 
the spins enter in a highly correlated quantum state, which in the case of classical spins 
is called a "cooperative paramagnet" 0, and a possible spin freezing below Tp into a 
Spin Glass phase [[7|, ^ ^; 2) a shift of spectral weight to low energies from high energies 
of the order of Qcw 3) Schottky-type anomalies in the specific heat in case there 
are no states in the singlet-triplet gap and a specific heat oc T" with singlet states in 
the singlet-triplet gap, while the susceptibility is always thermally activated ; 4) small 
external perturbations can lift (part) of the ground state degeneracy and thereby may 
radically change the type of (exotic) ground state: for example lattice distortions related 
to a Jahn-Teller effect |0 |12| bond disorder |T^, field induced order [1^, |1^ or the 



formation of a spin glass phase; 5) order-by-disorder: such small perturbations can also 
be quantum and/or thermal fluctuations (entropic ground states selection) |18|, pJ| ; 
6) anomalous, non-universal, critical properties 7) peculiar spin dynamics 8) in 
case of Ising spins quantized along the local < 111 > axis a cubic spin-ice type ground 



state p2| , |23| , which, in case that the dipole-dipole interactions are stronger than the 



exchange interactions, is described as Dipolar-Spin-Ice |2^, |2^, p7 | 

Reviews of experimental work on the susceptibility of Rare-Earth Stannate pyrochlores 
can be found in pH], for the magnetic properties of compounds with a pyrochlore spin 



lattice in |29|, and for the bulk magnetization of Rare-Earth Titanate pyrochlores in [3C, 



^. Reviews of theoretical work are given by Bramwell and Gingras Moessner 
and a review by Lhuillier and Misguich where the pyrochlore system is discussed in 
the perspective of the general class of frustrated quantum spin systems. 

Anderson [|33| started the study of order on the pyrochlore lattice and compared it 
with the positional disorder of H in hexagonal ice. He predicted the absence of Ising- 
type order down to T = in case of antiferromagnetic nearest-neighbour interactions for 



spins quantized along a tetragonal axis, calling such a system "Cubic-Ice" p^ , |35[| . We 
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can sum up the present understanding of the Classical Antiferromagnetic Heisenberg spin 
s = 1/2 on the pyrochlore lattice as follows: this spin system does not order from the point 
of view of classical spins: the highest J(q) modes are two-fold degenerate without any 



dispersion ^ meaning that below the mean-field ordering temperature Tc^mf = 
'^/'iJnn the spin system moves in these modes, thereby using only 1/2 of the total number 
of degrees of freedom. Moessner and Chalker ^ showed by an explicit construction 
procedure that all ground state configurations are connected, without any energy barrier. 
However certain correlations remain, which give rise to certain features in the neutron 
scattering: a few long-range correlations in high-symmetry directions is build into the 
ground state. 

Various authors have studied possible "order-by-disorder" in this lattice, introducing 
thermal or quantum fluctuations on top of the classical ground state, see for example |T^ . 
Most authors, probably starting with Reimers ||3^, first diagonalize the Heisenberg an- 
tiferromagnetic spin Hamiltonian with nearest-neighbour exchange on an isolated tetra- 
hedron and then construct the pyrochlore lattice with these non-overlapping tetrahedra, 
and introduce the interaction between these tetrahedra as a perturbation: XJnn- The 
tetrahedron-tetrahedron interaction has been calculated up to 3rd order in A. Garcia- 
Adeva and Huber ^Tj] used the most simple approach by neglecting these tetrahedron- 
tetrahedron interactions: in this way long-range order is excluded from the outset. The 
asymptotic "Curie- Weiss" temperature found is 1/2 of the value obtained from the canon- 
ical mean field approximation. Canals and Lacroix [42, 53| performed a density matrix 
expansion up to and found that at this order the degeneracy is weakly lifted (oc 10^^ J) 
giving rise to a coUinear AF spin structure: the spins along < 110 > are parallel while 
nearest-neighbour perpendicular < 110 > crossings are antiparallel. The correlation length 
is at most of the order of one lattice spacing. Their approach can not elucidate the char- 
acter of the ground state: whether it has long range order in one of the exotic ground 
states mentioned above. 

Several authors have discussed this possibility by calculating the tetrahedron-tetrahedron 
interactions in the restricted space of the two-fold degenerate singlet states and 3-fold 
degenerate triplet states of the tetrahedron. The latter are essential for inter-tetrahedron 
interactions. Although starting from different points of view, the ground state of this spin 
system seems to be some singlet -dimer covering of the pyrochlore lattice: Tsunetsugu Q 
used a chiral ansatz and found vanishing chirality: the ground state is a four sublattice 
system of tetrahedra with singlet dimer order on three sublattices each occupied by one 
of the possible tetrahedron singlet dimer states, and the 4th sublattice is locally disor- 
dered. The ground state of the singlet dimer lattice can be characterized as a Valence Bond 
Solid. Singlet excitations will appear in the singlet-triplet gap: these are oscillations in 
the dimer pattern. A similar ground state has been found by Eljahal et al [|^, by Isoda 
and Mori [Q. A dimer covering was already discussed by Harris et al |^^. Recently 



Berg et al pushed this approach one step further and used the two-fold degenerate 
singlet ground state of a super tetrahedron - consisting of 4 tetrahedra connected each by 
different bonds - as a basis for a perturbation calculation. To each super tetrahedron is 
assigned a pseudo spin, these super tetrahedra form a bcc lattice. From a Monte-Carlo 
study, treating these pseudo spins as classical entities. They find that the ground state is 
ordered anti ferromagnetically, with ferromagnetic layers. 
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Adeva-Garcia and HuberEDI have studied the magnetic properties of spins on the 



pyrochlore lattice using the constant-couphng approximation and recently applied it to 
ZnCr204P]. The same authors studied the critical properties of these pyrochlore spin 
systems using a Renormalization Group appraoch||5^ . 

When one looks into the pyrochlore structure one finds that each atom on this lattice 
shares six planar hexagons, two in each < 111 > direction, of which each of the three pairs 
of hexagons are edge sharing. The hexagons in a (111) plane form a Kagome lattice. So 
this pyrochlore lattice also may be characterized as a lattice of inter-penetrating Kagome 
lattices. 

Let us now turn our attention to nn ferromagnetic interactions. There is at first sight 
no reason why the spins on a lattice of GST would not order at finite temperature. First 
consider the case that the spins are Ising spins pointing along the local < 111 >. In that 
case there are six degenerate states with a two-spin-in two-spin-out configuration on each 
tetrahedron [^, and the spin states can be mapped on a model describing the two 



possible positions of Hydrogen in the ice structure [0]. This leads to so-called ice rules 
for the possible low temperature spin configurations. A clear exposition of these rules and 
its consequences can be found in P2|, Bl, E^. Because of this ground state degeneracy. 



the entropy at T = K is finite, and the macroscopic ground state of this spin system is 
called spin-ice. 

In the present paper we apply the Gluster Variation Approximation (GVA) for the con- 



figurational entropy, in which one divides the lattice in a set of overlapping clusters 1 52 



(basis cluster) to study the magnetic properties of Quantum Heisenberg spins on a 
pyrochlore lattice. We use the tetrahedron as basis cluster. We include tetrahedral 
anisotropy, and an anisotropy in the exchange interactions caused by a tetragonal dis- 
tortion of the lattice. Single-ion anisotropy and dipole-dipole interactions are discussed 
briefly. 

Recently this approach has been applied to the problem of classical Ising spins on 
a pyrochlore lattice by Yoshida et al The properties of this classical Ising model 

are quite different from the ones we report in this paper, especially the ground state 
degeneracy. 

The GVA is simply the Mean Field approximation in case the basic cluster is a single 
site. In case the whole crystal is taken as the basic cluster the approximation becomes 
exact. In between one has to make an intelligent choice for the basic cluster: we have 
chosen the tetrahedron as the basic cluster. The tetrahedron is the smallest unit possessing 
all the symmetry of the lattice, it is the unit which seems to contain the basic magnetic 
properties of the spin system: the high frustration against order in case of AF order. 
Another appealing property is that one can easily diagonalize the spin Hamiltonian for 
such a cluster, and the Glebsch-Gordan coefficients are also easily calculated. 

The hexagons next to the triangles, are the smallest self-avoiding closed loops on the 
pyrochlore lattice. So any description based only on the tetrahedral clusters, means we 
neglect these hexagon closed loops. In this sense our approximation is similar to that for 
a Husimi tree lattice, although our approximation does not suffer from the rather 
peculiar increase in the number of surface clusters, see also ||56| . 

This paper is organized as follows. In section |^ we introduce and discuss qualitatively 
a Hamiltonian describing the spins on this pyrochlore lattice. We transform the site 
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Hamiltonian to a Hamiltonian for tetrahedra ( |2.lD . In section we discuss the spectrum 
of a 16 site cluster. We then derive in section p.3| the susceptibihty for a pyrochlore lattice 
covered with non interacting tetrahedrons. In section p.4| we give some basic details of 
how we apply the CVA to this Heisenberg Quantum system. We derive an expression 
for the (linear) susceptibility ( p.4.2|) . In the next section ^ we present first the critical 
temperature (|3.1|) and we present the zero temperature phase diagram for ferromagnetic 



and AF nn interactions (|3.2|) , and then present the principal numerical results on the 
susceptibility, magnetization as a functioon of magnetic field and temperature ( |3.3| ). In a 
final section ^ we give a brief summary and conclusions. 

2 THEORY. 

The general Hamiltonian to describe spin-spin interactions is the Heisenberg Hamilto- 
nian. We consider nearest-neighbour (Jnn), and briefly we will also discuss qualitatively 
the consequences of next -nearest-neighbour (J„„n), and next-next-nearest-neighbour 
{Jnnnn) cxchaugc interactions. We assume that each of these interactions does not depend 
on the details of the interaction path in the pyrochlore structure: nn, nnn, and nnnn in- 
dicates by how many cation-cation links these cations are separated. From an analysis 
of ESR spectra of Cr^^-doped spinel ZnGa204 Henning derives the Cr^+-Cr^+ ex- 



change constants on this spinel B sublattice. These exchange constants show only small 
deviations - of about 10 % - from the average value of each type of exchange path. We 
also consider the possibility of local single-ion trigonal and tetragonal anisotropy. In case 
of the single-ion anisotropy the local trigonal field is along a local < 111 > axis pointing 
to the center of the tetrahedron. The tetragonal anisotropy can be due to a tetragonal 
distortion, due to uni-axial pressure, of the cubic pyrochlore or spinel structures. It will 
be accompanied by a change in the two nn exchange interaction in the {100} planes with 
respect to the other four nn exchange interactions in the {111} planes (see figure |^). The 
Hamiltonian for this spin system in a magnetic field H reads: 



H J nn ^ ^ ' ~^ J nnn ^ ^ ' ~^ J nnnn ^ ^ ' ~^ ^ ^ ^ ' 
i,j i,j i,j i,j;pairs 

+Dos E« ' - 34') + Dt E(S' - 35^ J + E gf^sH ■ s„ (2.1) 



The fourth contribution is a sum over the exchange interactions in the {100} planes. We 
reserve capitals Sq, for the size of the cluster (tetrahedron) spins and small s for the local 
atomic spins. The magnetic field H in the last term is along a tetragonal axis (see figure 

i- 

The s' in the single-ion contribution Dqs term indicates quantization of the spin along 
the local trigonal < 111 > axis. In case the single-ion anisotropy Dqs > the local < 111 > 
axis is the easy axis, while in case Dq^ < 0, the spin rotates in the x'y' plane perpendicular 
to this axis. In case Dq^ > and much larger than the exchange interactions it is often 



argued that one may approximate the spins as classical (Ising) spins |5^. In case of a 



tetragonal distortion along one of the 5*4 axes of the tetrahedron one finds 

H,, = iDo.(3cos2^4-l)E(«'-34) (2.2) 
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cos ^4 sin ^^4 + Si^Siz — S2zS2+ — S2+S2Z + i{s3zS3- + S^^S^z — S4ZS4- — ■54_S4^ 



3 

-Dqs sin^ 6*4 



[si - si) + {si - si) - {si - si) - {si - si: 



where 64 is indicated in figure 0. The single-ion anisotropy vanishes in case s = 1/2. 

Next let us briefly discuss the tetragonal anisotropy Dt- When the tetrahedral 
anisotropy Dt > 0, the tetrahedron spin S is aligned along a tetrahedral axis, and in 
case Dt < 0, this spin is rotating in the plane perpendicular to a tetrahedral axis. This 
term can also be written in terms of the single-site spins: 

{Jnn + Dt) Y1 SiSj - 3Dt Y1 SizSjz + '^Dt "^{s^i - 3sl = 

<i,j> <i,j> i 

+ 2DtY.{sI-3sI). 

<i,j> <i,j> i 

(2.3) 

The tetrahedral anisotropy interpolates between a Ising model (J„„ = —Dt) and an 
XX Z model to a Heisenberg {Dt = 0) and XX model {Jnn = '^Dt), by changing the 
ratio DT/Jnn- In the next section we will get more insight in this term by writing the 
whole Hamiltonian in terms of tetrahedron spins. 

2.1 The Tetrahedron Approach. 
2.1.1 The Hamiltonian. 

For the CVA with the tetrahedron as a basic cluster we need the spin states of this 
tetrahedron. We define for tetrahedron a the total spin as = Si + S2 + S3 + S4. For the 
A term we introduce L = Si + S2 and R = S3 + S4. One can now rewrite the exchange 
part and magnetic field of the original Hamiltonian in terms of these tetrahedron spins 

H {Jnn '^Jnnn ~l~ 3Jnnnn) ^ ^ S^, -|- {Jnnn '^Jnnnn) ^ ^ Sq. ■ -|- Jnnnn ^ ' ' 

nnn 



+ W ^/^bH.S. + a ^ (l^ + Kl) + Dt E(S^ - SSl) + H,, 

^ a a a 

(2.4) 



The sum is over all tetrahedra in the lattice, including the ones which share a corner: 
so in case of the pyrochlore lattice there are N/2 GST. The sum over a and (3 is over 
all tetrahedra under the condition that they are nn or nnn tetrahedra. We have dropped 
terms of type JY^i sf , as these only add a constant energy to the Hamiltonian. The first 
term gives the energy levels of an isolated tetrahedron, with an effective exchange in one 
tetrahedron of J„„ — 2J„„„ + 3J„„„„. This correction on J„„ should only be taken into 
account when also the the corresponding nn and nnn are taken into account: the second 
and third terms give these interactions between the total spins on nn and nnn tetrahedra. 

We note that several of these contributions {Dt, A) also appear in a summation of 
the dipole interaction within a tetrahedron: 



D 



D 



D 



3s: 



2D 



D 



+ - 3{Ll + Rl) +■■■ (2.5) 
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where D^, = g'^fi'^/R^, R is the interatomic distance in a tetrahedron. The dipole inter- 
action, diagonal in the tetrahedron cluster states, is taken into account (except a contri- 
bution of the form Ll + Rl), the last term in eq. p.5| . All single and double spin flip terms 
appearing in the dipole-dipole interaction are not taken into account. These contributions 
cause the dipole-dipole interaction neither to commute with Ms, nor with S"^. 



2.2 The eigenstates of five coupled tetrahedra. 

In order to discuss the ground state of a lattice of CST, we calculated the ground state of 
a finite cluster consisting of 4 outer tetrahedra each sharing a corner with the same central 
tetrahedron: 5T super-tetrahedron. First let us briefly consider the energy levels of the 
tetrahedron. These are in figure |^ and table In case of nn AF xchange interactions 
the ground state is a degenerate pair of singlets and in case of F interactions a quintet 
state S = 2. This pair of singlet states is split by the + R^ perturbation into two 
singelts, one with L = 0, R = 0: the two dimer singlet state, and the other L = 1, R = 1: 
the plaquette or tetrahedron Singlet state. Also the dipole-dipole interaction splits this 
singlet. 

We calculate the ground state energy in case of AF interactions. We proceed as fol- 
lows: the Heisenberg Hamiltonian is rewritten in terms of the tetrahedron spins Si, for the 
outer 4 tetrahedra. The eigenstates can easily be found using the Clebsch-Gordan cou- 
pling coefficients. These states are represented by the following set of quantum numbers: 
\Q, Mq, Di, D2, Si, S2, S3, Si, Li, Ri, L2, R2, L3, i?3, L4, Ri), where Q and Mq are the total 
spin and magnetic moment of the outer 4 tetrahedra, Di is the spin of the tetrahedra 
with Si and 5*2, and D2 of 5*3 and S^. The Li and Ri are defined above. This basis is di- 
agonal in the Hamiltonian for the exchange interactions of the outer four tetrahedra. The 
6 exchange interactions in the central tetrahedron have now to be considered separately. 
For Mq = 0, and Q = 8, 7, 6, 5, this can be done easily. For the other Mq values we 
have chosen the following scheme in order to determine the ground state and first excited 
states: we choose a maximum value for Si, Di and Q, for Mq = and restrict the sum 
over all tetrahedron spins Si. That is, we restrict the number of singlet-triplet-quintet 
excitations in the cluster. Our results for the ground state energy Egs and singlet-triplet 
gap AsT are in table ^ We find that the ground state is a set of 16 singlets. There is 
only coupling between the four outer tetrahedra by singlet-triplet excitations. We find 
Egs = —12. 82 J for our T5 super-tetrahedron, which is — 0.81J/site. This value is smaller 
then obtained for the 16 site super-tetrahedron of Berg etal [Q. The singlet-triplet gap 
is about AsT = 1-56 J, which means a decrease of about 25 % with respect to indepen- 
dent tetrahedra. In table ^ we have collected some data from recent calculations of these 
quantities. 

The ground state degeneracy of the T5 cluster is determined by the 2-fold degenerate 
ground-state of singlets of each of the 4 outer tetrahedra. The state on the central 
tetrahedron is completely determined by the state on each of the four outer tetrahedra. 
When we extend this description to the complete pyrochlore lattice - which is from the 
point of view of the tetrahedra a diamond lattice - the ground state is 2^/^-fold degenerate 
for a lattice of N sites in case of AF interactions. 

In our calculations on the T5 cluster we observed a partial lifting of the degeneracy 



7 



for certain choices of the restriction on the maximum number of excitations (the size 
of Hilbert space): to get the correct degeneracy one has to diagonahze the complete 
matrix. In this T5 cluster we do not find lifting of the degeneracy due to inter-tetrahedron 
interactions . Diagonalization of only part of Hilbert space gives incorrect degeneracies 
and spurious eigenvalues. 

Such singlet gap-states seem to appear as a function of the cluster size, or in a real 
system as a function of the spin-spin correlation length. This implies that the extent of 
spin singlet states depends on temperature, and the singlet-singlet excitation spectrum 
collapses to a single highly degenerate level or a narrow band (-C J) of singlet states: 
that is, transfer of spectral weight from high to low frequencies with decreasing tempera- 
ture. The appropriate correlation lenght if this quantum order increases with decreasing 
temperature. At high temperature (T > J) this correlation lenght is of the order of one 
interatomic distance. 



2.3 The paramagnetic susceptibility of independent tetrahe- 
dra. 

With nn exchange interactions J only, and neglecting the restrictions due to corner sharing 
of spins, one derives easily the following expression for the susceptibility PH] : 

_ j T.s 9sS{S + 1){2S + 1) eM-Es/kBT) 
^ ksT Es9s{2S + l)exp{-Es/kBT ' 

where gs is the degeneracy of a tetrahedral S'-state excluding the {2S + 1) degeneracy 
which is included explicitly in this expression. The energy of each S state is given by 

Es = JSiS + l). (2.7) 

In the high temperature limit one derives: 

^""^ ^ 3kBTl-{Gi-G2)J/kBT' ^^-^^ 



where 



Es9s{2S + l){S + l)S 
J:s9s{2S + 1) ' 



^ j:s9si2S + l){S + lfS' 
J:s9s{2S+1){S + 1)S ■ 



One can now easily calculate the asymptotic Curie temperature: 

kB^cw _ 1 
J Gi — G2 



(2.10) 



This asymptotic Curie temperature is not related to any long range order. For the case 
of s = 3/2 - for example Cr^"*" with d^ - one finds an asymptotic Curie temperature given 
by: 

kBOcw = -^J, (2.11) 
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which is 1/2 the value one would obtain from the usual Mean Field approximation! One 
finds the same factor 1/2 for the Kagome lattice in case one treats the lattice as consisting 
of independent corner-sharing triangles. 



2.4 The Cluster Variation Approximation with tetrahedron as 
basis cluster. 

In this section we give some details of our calculation of the probability functions for 
quantum spins for the case of a tetrahedron cluster with s = 1/2. The calculation for 
higher spin values is only slightly different. The CVA is often applied to the classical 
problems of order-disorder transitions in alloys, and the problem of Ising spins. In the 
present case we apply it to the problem of Quantum Heisenberg spins. The approximation 
we make is that there is no coherence between the quantum states on different tetrahedra. 
At high temperature this leads to the classical Mean Field limit. However we will see that 
in certain cases below a certain temperature T5 the CVA breaks down: the entropy in 
this approximation becomes negative for T < T^. 

The states on the tetrahedron are written as linear combinations of the determinants: 
\mim2Tn^mi\, where is the spin moment on i. The spin eigenfunctions can be written 



as: 



|\E'(S', M, L, i?)) = |\I'(S', M, a)) = ^ (mi, m2, ms, m4|S', M, a)|mi, m2, ms, 7714), 



(2.12) 

where (mi, m2, ms, m4|S', M, a) are the Clebsch-Gordan coefficients. The a's indicate 
one of the eigenfunctions of the possible degenerate set of eigenstates with spin S and 
momentum M with different pairs {L,R). The restrictions are \L — R\ < S < L + R, 
M = Mr + Ml, < L < Si + S2, < < S3 + S4, Ml = mi + m2, Mr = nis + m^. 
The probability that a certain site has magnetic spin momentum m is: 

(2.13) 

S,M,oi 

where ys,M,a is the probability that the tetrahedron is in one its eigenstates \S,M,a), 
and am-s,M,a is the sum of the square of the (Clebsch-Gordan) coefficients of the states 
|mi, m2, ms, m^), with mi = m in the expansion of the spin eigenstates of the cluster in 
terms of these determinants |2.12| . For degenerate states it is important to symmetrize 
these coefficients, giving equal weight to each of these states. 
The total average magnetic moment per site is now: 



(m 



m 



We define: 

bm;S,M,a = 0'm;S,M,a — flm;0,0,«o5 (2-15) 

where the index ao is (L, R) = (0, 0). These 6-coefficients are in table ^. We define: 

(2.16) 
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The free energy of a system of A^^ clusters on a lattice can be written: 



F = NtY.[ij,E^-kBT\TiW{^)i 



(2.17) 
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where is the energy of a cluster in state 7 = (S", M, a), and Wi^-^) is the number of 
ways to distribute a set of states {7}, each occurring with probability c<j^ over the lattice. 

In case of the tetrahedron as basis cluster, the number of ways to distribute these 
tetrahedron states over the lattice of CST is calculated as follows: first distribute these 
tetrahedron states over the whole lattice. Zt tetrahedra will share the same corner. So 
there are Nt = Zt/AN tetrahedra, each characterized by a spin eigenstate on the lattice. 
The number of ways of putting these tetrahedron spin states on the lattice is: 



In order that a certain distribution belongs to a realizable one, the local spin moment 
for each of these tetrahedra should be the same on the shared atoms, called the overlap 
cluster. In case of the pyrochlore lattice the overlap cluster is the point cluster, and 



This restriction on the distribution of these tetrahedron spin states can not be applied 
rigorously, but on the average: the CVA for the entropy. At high temperatures the spins 
are distributed over all spin eigenstates of the tetrahedron: this is the classical regime of 
spin fluctuations. At very low temperature, when the tetrahedron spins are distributed 
over a small number of tetrahedron spin states, one is in the regime of quantum spin fluc- 
tuations. In case no classical long-range spin order interferes this classical approximation 
of the entropy can break down. In case the degeneracy of the tetrahedron ground-state 
is too low, the entropy within the CVA is negative below a certain temperature. When, 
at low temperature, the tetrahedron ground-state is a single non degenerate singlet state, 
no correction is necessary: due to quantum spin fluctuations the average spin on each site 
vanishes: there is self-averaging. The description of such a state is outside the validity 
regime of the CVA. We take into account correlated spin fluctuations of the size of the 
tetrahedron, correlated spin fluctuations over larger distance are not described within the 
present approach. With this self-averaging one might argue that in a certain temperature 
range before the system enters the regime of correlated spin fluctuations, the state on nnn 
tetrahedra are independent and one could approach the susceptibility by eq. If that 
is the case the asymptotic Curie-Weiss temperature found from magnetic susceptibility 
data at relative low temperature gives an exchange interaction 1/2 the actual value. 

The correction within the CVA is now calculated as follows. Distributing Nt tetrahedra 
over the pyrochlore lattice generates ZtN spins. On each site of the pyrochlore lattice 
there are Zf = 2 spins, which we call pseudo sites. First calculate the number of ways 
of putting ZfN spins with at any moment either spin up or down on the lattice of ZtN 
pseudo sites. Of these arrangements, only a fraction is correct: all spin momenta on the 
Zt pseudo sites belonging to the same site are the same. Now the average fraction of 
correct point configurations is: 



Q 



Nt\ 




U-y{Nty^y. 



Zt 



2. 



1 l-Zt 



r 



(2.19) 
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Combining this correction with the number of arrangements of the tetrahedron configura- 
tions on the pyrochlore lattice we obtain for the average number of correct arrangements 



1 


Zt/4 


1 


l-Zt 






_(a;_i/2iV)!(a;i/2iV)!_ 





(2.20) 



In case of s > 1/2 one should replace the product (a;_i/2-/V)!(a;i/2-/V)! by the appropriate 
product for spin s: l\m=-s,si^mNy.. 

Substitution of this equation for the total number of spin configurations in the expres- 
sion for the free energy and using Stirling's formula for large we obtain for the free 
energy per site 

f = F/N = Zt/Aj2y,E^ 



Zt/4:^y^ln{y^) + {1- Zt) ^ x„ln(x^) 



(2.21) 

In order to obtain the equilibrium values for the probability functions we have to 
minimize the free energy with respect to these functions y^, and find: 







dl 
dy^ 



Zt/4. 



[E^-Eo)+kBT\\n^ + ^ 

yo 



E 



dXr, 

dy^ 



ln{xm) 



(2.22) 

where we used the fact that these distribution functions x and y are normalized, yo and 
Eq indicates the state with S — 0, L — 0, R — 0. Using the definitions of the coefficients 



we can write: 

dy-, 



We eliminate the state from these equations and write for 



the energy of the eigenstates in a magnetic field: 

Es,M,a = E%^M,a + gfJ'sMH. 

where Eg j^^^ — Eg _j^^^. The equilibrium equations become: 

^^iMi^ = -Qa.2C,B^,s.M.. exp(2MVT), 

yS,-M,a m 



(2.23) 

(2.24) 



where we defined: h — g/iBH, and we set ks — 1 and the constant Ct = 2{Zt — 1)/Zt. 

In case the spin system is non-magnetic - no long-range order - one can easily solve 
these equations. The distribution functions are independent of the magnetic quantum 
number M, and the probability that a spin is in one of its 2s -|- 1 magnetic moment states 
m is in case of Heisenberg spins: l/(2s -|- 1). 

In case of finite single-ion anisotropy, the Xm should be evaluated using the appropriate 
probability distributions for the non-magnetic solution. 

The distribution over the tetrahedron states is in the non-magnetic case without 
magnetic field given by: 

ys,M,a = Z/oexp(-(£'s,M,a " Eo)/kBT), (2.25) 

and 

yo = 1 - J2 ysM,a- (2-26) 

S,Ma^O 
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2.4.1 The entropy <S for s = 1/2. 

The free energy is for s = 1/2 in the disorder regime: 

F ^ ~kBT In Zo + keT ln2 (2.27) 
and the entropy in this regime is: 

5/A;« = lln(V4) + ^r^^ (2.28) 

where 

Zo = ^exp(-£;^/A;Br) (2.29) 

7 

In the disorder regime the single-site spin momentum probabihties x„i=±\i2 = 1/2 are 
constant. For T ^ the last contribution to the entropy can be neglected. We find that 
for < 4 the entropy becomes negative: the restrictions imposed by the CVA are too 
strong. The temperature where this occurs we call T5. We find that in case there is no 
long-range order at any finite temperature, then in case the tetrahedron has a ground 
state with a degeneracy less than 4, the spin system has in some way to enter a state of 
correlated spin motion in order to satisfy the thermodynamic requirement that 5 > 0. 
Already far above Tg one should observe strongly correlated spin fluctuations extending 
over more than one tetrahedron. 

It seems that the restrictions to determine the number of realizable configurations 
imposed by the CVA are too strong. At low temperatures, when only singlet tetrahedron 
states are occupied, there is self-averaging and no CVA-type restrictions are necessary. 

2.4.2 The magnetic susceptibility. 

In order to calculate the magnetic susceptibility we linearize the functions which describe 
the distribution of the spin configurations in the magnetic field H and the magnetization 
M. as follows. First we write: 

yS,±M,a = yO;5,M,a(l ± rns,M,a), (2.30) 

where ms,M,a is a new function linear in the magnetization and magnetic field. In the 
following the summations over M arc restricted to M > 0. We linearize the exponent 
which contains the magnetic field and get: 

{l + 2ms,M,a) = il + 2Mh) n (2.31) 

m=—s,s 

Next we write an expression for the single-site distributions in terms of these cluster 
distributions, and we find 

= ^0,s,m + XI Bjn.s,M,ayS,M,atTT'S,M,a, (2.32) 
S,M>0,ajtO 
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where we defined 

Xo,s,m— 51 bm-S,M,ayS,M,a + bmiO)- (2.33) 

We have now reduced the number of unknown state distribution functions to the states 
without taking into account the M degeneracy. For s = 1/2 there are 6. 
One can easily see that Bm = —B^rn- So we can write {h — H/kBT): 



l + 2m5,M,as = (1 + 2M/1) 

■■r-r r D ^^ "1 Bm-S,M,a 

llm>0 ^Os + Z^5,M>0,a^0 Bjn-S,M,ayS,M,a'fnS,M,as) 



Ylm>Q ^Os ~ S5,M>0,a^0 Bm;S,M,ayS,M,a'>TT'S,M,a) 



(2.34) 



Tliis can now easily be evaluated to lowest order in the magnetic field and the magneti- 
zation for each m state. We define the following quantity: 



m 



m 



XI Bjn-S,M,ayS,M,a'^S,M,Oii (2.35) 



S,M>0,ai^Q 

and we can now solve for m^: 



2 2 

— K^h — Ct Yj^nim — Ct X^ YmmilT^mn (2.36) 

where we defined the following functions: 

Km — 53 Bjn;S,M,aMys,M,a, (2.37) 
S^O,M>0,a 

Bm2;S,M,ayS,M,a- (2.38) 

S7^0,M>0,a 

In the cases s = 1/2 and s = 1 this set of equations can be solved easily. In the case 
s = 1/2 one finds 

mi/2 = i^i/2/(l + Ct^—Y,/2)h, (2.39) 

3^0,1/2 

with = 1/2. 

The magnetization can be written as: 

M = J2 = 2 53 mmm (2.40) 

m m 

The second equation for A4 follows easily from the definition of and m^. The factor 
2 appears, because the sum in over M is only for M > 0. 
The magnetic susceptibility is now given by: 

^ = M/H^'^^^^^^^^ (2.41) 
H 



and in case of spin s — 1/2 we get 



X = ^ = , , (2.42) 
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The functions K and Y can easily be calculated from the Clebsch-Gordan coefficients 
and the probability distribution functions for the non-magnetic system. We find that the 
asymptotic Curie- Weiss temperature in the CVA is equal to the MF result. 

For s = 1/2 we find for the susceptibility in case of zero uniaxial fields {Dt = A = 0): 

X = -f 2 + 6e2^/^ • ^^-^^^ 

where we have set fc^ = 1. A comparison of the CVA result and an exact high-temperature 
expansion shows that the CVA 



X 



^ ^1 + 3x + + 



T 

already deviates at second order in x = J/T 



X 



1 + 3a; + 9a;^ + ■ 



^2.44) 



(2.45) 



T 

We are not aware of any high-temperature expansion for the CST lattice, only for the 



Kagome [31 



The magnetization is given by the solution of the following implicit equation in J\4.: 



(Ai 2 + i^i 1 )e T + ^ e^^^ e^T + ie^T 

M = ■ = — — = - (2 46) 

. 4A , . , . 4_ 6(J+0)+4A , o J+Pt nA .A\ l^'^^/' 

1 + + {AI2 + At^ + l)e T + (^+ ^ + 1)6^^?^ [2e^T + e^T 



where we defined: x = (1 + M)/{1 - M), Af^ = (xe^^/^ ± le-'^H/T^^-uoT/T ^ 
A± 2 = (Vie^/^ ± -^e-^/^)e-3^-/^. 

3 RESULTS AND DISCUSSION. 

3.1 The critical temperature for Heisenberg spins. 

We determine the critical temperature as the temperature where the susceptibility di- 
verges. The results are in table |^. In case s = 1/2 and only Heisenberg interactions, 
there is no magnetic long-range order at any finite temperature, neither for F nor for AF 
interactions. In our CVA to the pyrochlore lattice, AF Heisenberg spin interactions never 
give classical long-range order. In case of ferromagnetic Heisenberg interactions the CVA 
values are about 60 % of the usual Mean Field values, slightly increasing - as expected - 
for larger spin. 

We calculated for s = 3/2 the critical temperature for F Heisenberg exchange from 



a high-temperature approximation |62] with 6 coefficients using the Fade extrapolation 
and find Tc/J = —7.892. Compared with this, in principle exact calculation, the CVA 
critical temperature is still too high, indicating that spin fiuctuations are stronger than 
accounted for in our application of the CVA. Note also that our critical temperatures 
for ferromagnetic Heisenberg interactions are appreciably lower than those for the same 
coordination number using the interpolation equation of Rushebrook and Wood 
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We studied the effect of increasing the number of tetrahedra sharing the same corner. 
We found that any increase however small will lead to a finite critical temperature for 
s = 1/2 in case of ferromagnetic Heisenberg exchange. This is illustrated in figure |^. 
A similar conclusion was reached by Pinettes et al who - using MC simulations - 
analyzed the ordering temperature going from a pyrochlore to a fee structure. They 
found that for very small values of the ratio of the extra exchange interaction - needed 
to transform the pyrochlore lattice into the fee lattice - and the exchange interaction in 
the pure pyrochlore lattice, the system shows a coUinear long-range order. 

3.2 Ground state Phase Diagram for s = 1/2. 

In order to derive the ground state phase diagram for s = 1/2 we calculated the energy 
levels of the spin states of the tetrahedron, which can be found in table |l] and are illustrated 
in figure ^. The phase diagrams in the drS plane, where dr = Dt/J and St = ^t/J, 
are in figure | for ferromagnetic as well as antiferromagnetic interactions. We note that 
for various values of the distortion the ground state is non-magnetic. The lines separating 
the various regions can easily be found from the expressions for the energy. 

In case of F interactions we find the following equations for the boundaries between 
the various ground states: 

= 0, M = 0)-{S = 2, M = 0): 25 + Sdr + 3 = 0; 
(5 = 0, M = 0)-ls = 2, M = 2): 25 - 3c/t + 3 = 0; 
ls = 2,M = 2)-ls = 2, M = 0): dr = 0. 

So we find 3 regions: one region where the spins show long-range order (S* = 2, M = 2), 
and two regions where we find a non magnetic ground state: one where the ground state 
of the spin system is build from singlet dimers [S = 0,M = 0,L = 0,R = 0), and one 
build from the nonmagnetic M = component of the quintet state, indicated by RVT2. 

In case of AF interactions we find the following boundary equations: 
(S=2,M=2)-(S=0,M=0): 26 ~ Sdr + 3 = 0; 
(S=0,M=0)-(S=2,M=0): 26 + Sdr + 3 = 0; 
(S=0,M=0)-(S=2,M=0): dr = -1; 
(S=2,M=2)-(S=0,M=0): dr = 1. 

In this case we find four regions: the region with a classical long-range AF order, 
and three regions without magnetic long-range order, the region indicated by RVB in the 
diagram, with a ground state build up from singlet dimers, a region characterized by a 
singlet tetrahedron ground state with compensating triplet spins L = 1, R = 1, indicated 
by a Resonance Valence Tetrahedron state: RVTl, and a region where the ground state 
is build from the nonmagnetic spin moment M = of the tetrahedron quintet state, 
indicated by RVT2. 

In case of AF interactions and dr > 0, the ground-state becomes magnetic for dr > 
dT,c = (4s + 1)/ (8s + 1). In case of s = 1/2 the critical value dT,c = 1 and for very large s, 
dT,c 1/2. If we translate this to the case that dipole-dipole interactions are the source of 
the tetragonal anisotropy, then = -Dd/8, and in the large s limit we find for the critical 
ratio of the nn exchange and the dipole-dipole interaction: J/Dd = 1/4. For smaller 
ratios the 5 = 2, M = 2 state drops below the singlet 5* = 0, M = state. Whether 
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this can be applied to the Ising systems discussed in a recent review by Bramwell and 
Gingras , is not clear. In case we make the necessary corrections for J; their exchange 
interaction is -1.5 times our exchange constant, and we find from their work the following 
condition for magnetic order: J/D <1 , which is of the same order of magnitude as we 
have derived. 

The boundary separating the magnetic from the non-magnetic solutions is not a sharp 
line: these two regions are separated by a re-entrant region where there are two critical 
temperatures. The extent of this region increases with increasing tetrahedron distortion 
like illustrated in figure |. In figure ^ we give the critical temperature ksTc/ J as a 
function of the tetrahedral anisotropy dx, for various values of the exchange anisotropy 
S. This re-entrant behaviour is due to the occupation of a magnetic state with increasing 
temperature i.e. an entropic effect. It is a type order-by-disorder [O transition: an 
increase in thermal spin fluctuations induces magnetic order. 

Koga and Kawakami |^ have also derived a phase diagram for the Quantum Heisen- 
berg pyrochlore system with AF interactions. They studied the influence of the exchange 
anisotropy A. When comparing their result with ours we might identify the plaque- 
tte state with our 5' = 0,M = 0,L = 1,_R = 1 state and their dimer state with our 
S = 0,M = 0,L = 0,R = state. For 6 = 0, for an undistorted tetrahedron these two 
states are degenerate. Whether the system with A = — J and AF interactions, i.e., the 
total L and R exchange interactions vanish, should show magnetic order is questionable. 
It is a very open 3D network with a very low coordination number of 4. 

In flgure |^ we compare the ferromagnetic order temperature T^/J as a function of 
tetragonal anisotropy (It = Dt/J for various values of the exchange anisotropy 6 = A/J. 
We find that Ts < Tc in case of finite Tc. 

In the regions of the phase diagram where the ground state has a degeneracy smaller 
than 4, there is a finite temperature Ts where the entropy in the CVA becomes negative. 
In figure ^ we present W -the number of states accessible by the spin system so that the 
entropy per site is S /ks = \\iaW - for the case of ferromagnetic exchange J = — 1 as 
a function of T/| J|, without any perturbations. Because the spin system does not order 
for T > 0, the entropy remains finite for all T > 0. At high temperatures the entropy 
approaches the value for a paramagnetic spin system (In 2). 

We surmise that below Ts the spin system enters in a state with strongly correlated 
quantum spin fiuctuations. The latter takes over from the slowed down thermal fiuctua- 
tions. This occurs in all regions without magnetic moment degeneracy (M = 0). 

Such a ground state degeneracy can also be lowered by a magnetic field. In figure ^ we 
present some results of Ts as a function of the tetragonal anisotropy for 5 = and 6 = —3. 
In the first case the magnetic field favours the state with extended spin fiuctuations for 
all dr's while in case of 6 = —3 these fiuctuations are suppressed with increasing field for 
small tetragonal anisotropy, while for large dx the quantum correlated spin fiuctuations 
are favoured by a magnetic field. 

We have determined the extend of the re-entrant region for s = 1/2 to 7/2. In fig. IIU 



we present the re-entrant region for 6 = for s = 3/2. In table ^ we give T^ax, Drp^max 
values which characterize the extend of the re-entrant region for s = ltos = 7/2: this 
re-entrant region increases with increasing spin. We expect to find similar behaviour in 
case of trigonal single-ion anisotropy. 
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3.3 Magnetization as a function of field and temperature: 
the F case. 



In figure |ri| we present the magnetization as a function of the magnetic field for ferro- 
magnetic exchange (J = —1) and s = 1/2 along the lines A-C [6 = —3 and B-D 5 = 1 in 
figure for T = 0.01 |J|, 1| J|, 5| J|. We can distinguish three types of behaviour, which is 
directly related to the ground state in the phase diagram. 

1. In the region of (ferromagnetic) order {S = 2, M = 2) the magnetization behaves in 
the usual way as a function of the magnetic field. We note that for the unperturbed 
case {6 = 0, dx = 0) only a very small magnetic field is necessary to magnetize the 
system. 

2. In the region characterized by(S' = 0, M = 0, L = 0, -R = 0), that is the region 
of the phase diagram where due to the 6 and dx anisotropy, the magnetic states 
are pushed far above the non magnetic ground state. Here we observe, at very low 
temperatures as a function of the magnetic field, a single first -order phase transition 
from a non-magnetic state to a completely ordered state: M. = J^sat for a field 
he = —3 — 26 + 3dT- This transition smears out with increasing temperatures. At 
T = this is what is called a Quantum Phase Transition. 

3. The region characterized by(S' = 2,M = 0,L = l,i?=l). In this region the 
magnetization as a function of the field shows two first-order phase transitions as 
a function of the field at very low temperature: one at = Sdj- and the next at 
hc2 = ^dx (Quantum Phase Transitions). The plateau in the magnetization is at 
Ai = l/2M.sat- Like in the previous case these transitions are smeared out with 
increasing temperature. 

Note that in various cases the magnetization increases nearly linear with the field, and in 
some case it changes from a convex to a concave function with increasing field. 

These magnetic field dependences lead to a rather peculiar temperature dependence of 
the magnetization. In figure |TT] we present the temperature dependence of the magnetiza- 
tion along the B-D line (5 = 1) in the ferromagnetic phase diagram, for positive dr values 
(S" = 2, M = 0, L = 1, i? = 1). Fundamentally we observe three types of behaviour of the 
magnetization as a function of the temperature. For small field the magnetization has a 
maximum as a function of the temperature for all values of > 0. With increasing field 
for small values of dx the system is ordered at low temperatures, while for large values of 
dx the magnetization has a maximum. For intermediate values of dx (see figure |TT]) the 
magnetization as a function of temperature exhibits a relative minimum. Note that in 
case this relative minimum is present in the magnetization as a function of temperature, 
than the zero temperature magnetization: A^(T ^ 0) ^ O.SA^sat, independent of the ap- 
plied magnetic field. This is related with the plateau M. = l/2Aisat in the magnetization 
as a function of field. 

In the region (S" = 0, M = 0, L = 0, i? = 0) one only finds one peak in the 
magnetization function of temperature. 

In figure |l^ we present the magnetization in a small field for 6 = —3 for various values 
of the tetragonal anisotropy dx in the re-entrant region. The peak in the magnetization 
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is clearly visible. Note that even outside the ordered region there is a finite magnetization 
which does not decrease to zero for vanishing field. 

In figure |n| we illustrate the behaviour of the specific heat as a function of temerpature 
in the re-entrant region for ferromagnetic exchange. Note that there appear two peaks on 
top of a Schottky peak. The behaviour of the specific heat at low temperatures is uncertain 
as the CVA breaks down (negative entropy). We conjecture that trigonal perturbations 
can induce similar beaviour. 



4 SUMMARY AND CONCLUSION 



We conclude from our calculations that in case of AF nn exchange only there is no 
classical long-range spin order on the pyrochlore or B sublattice of the spinel structure 
at any finite temperature. This is a consequence of the high degree of frustration on this 
lattice for AF interactions and its low coodination number. In case of the simple cubic 
lattice, which has the same coordination number (6) as the pyrochlore lattice, there is 
spin order below a certain finite temperature with nn interactions only. 

In case of F nn exchange only there is no spin order for finite temperatures for s = 1/2. 
For s > 1/2 one finds a finite critical temperature which is appreciably smaller than that 
found in the MF approach. However from a HTE using the coefficients in 
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one finds 

an even smaller ordering temperature, which indicates that even in the case of F nn 
Heisenberg-type interactions, spin order is frustrated. Such a frustration of spin order in 



case of F interactions has been discussed by Bramwell and Harris p2|, ^ for the case of 
classical Ising spins, each pointing along the local < 111 > axis of a tetrahedron: the 2 
spin in 2 spin out ground state configuration is sixfold degenerate, leading to a finite zero 
temperature entropy. In case of quantum spins there is no such degeneracy. The reason 
why the s = 1/2 pyrochlore system with F interactions does not order, is due to the low 
coordination number (6). One could also argue that from the point of view of the S = 2 
tetrahedron spin, the coordination number is only 4, and the smallest closed loop is a 
hexagon. Heisenberg spins on such a four-fold coordinated structure will probably not 
order. 

We did not include the effects of long-range dipole-dipole interactions, which recently 
have been proposed to be responsible for the so-called dipolar spin-ice ground state. 
Our appraoch is valid in case Jnn > Dd p.5| , so for systesm containing magnetic transi- 
tion metal ions like spinel with magnetic ions on the B sublattice ZnCr204, CsNiFeFg, 
CsMnFeFeU 



CsNiCrF 



6; 



GeCuoOJ|67f 



|68 |, than to systems with rare-earth 



_ LiV204 

iopns on a pyrochlore lattice. 

We have found in various regions of parameter space re-entrant magnetic behaviour 
when lowering the temperature: non-magnetic (Paramagnetic) —>■ Ferromagnetic — > non- 
magnetic (Two-Dimer-Singlet or Tetrahedron Singlet). The high temperature phase is 
paramagnetic, while the low temperature phase is some type of spin fiuctuating state. 
Such re-entrant behaviour has been reported as a function of the magnetic field for the 
garnet Gd3Ga50i2 by |Q, H, |70| and recently been discussed by Tsui et al|l^ . 

The magnetization as a function of the field depends strongly on the model param- 
eters. We find magnetization curves which show a single step for large field or a two 
step behaviour. Stepwise behaviour has been found for classical Heisenberg magnets by 
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Zhitomirsky et al |]7T|. See also the discussion of magnetization plateaus by Lhuillier and 
Misguicli [0. Such steps have been observed in various pyochlore systems p9|, [30[] . 



We have also calculated the specific heat for s = 1/2. Details of the calculation of the 
specific heat will be published elsewhere. For the zero-field specific heat of non-magnetic 
systems we find usually a single peaked Schottky-like anomaly, and in case of the re- 
entrant phase we find two peaks in the specific heat, possibly on top of a Schottky-like 
peak. Such peaks have recently been observed in the specific heat of Gd2Ti207 [|72| . 



Summarizing we find indications from this CVA of a rich magnetic behaviour for spins 
on the highly frustrated pyrochlore lattice as a function of tetragonal distortions: Two 
dimer singlet or tetrahedron ground states, re-entrant behaviour due to entropic ground 
state selection, magnetic field induced order Quantum Phase Transitions), together with 
steps in the magnetization as a function of field, Schottky like anomalies in the specific 
heat. 
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Table 1: The energy levels of the spins on a tetragonally distorted tetrahedron for s = 1/2. 



S M L R 


E(S,M,L,R) 


S M L R 


E(S,M,L,R) 








2 ±1 1 1 


6 J - 3D + 4A 


10 10 


2{J + D) + 2A 


2 11 


6(7 + 1?) +4A 


10 1 


2{J + D) + 2A 


11 


4A 


1 ±1 1 


2J- D + 2A 


1 ±1 1 1 


J-D + 4A 


1 ±1 1 


2J-D + 2A 


10 11 


2(J + D) +4A 


2 ±2 1 1 


6J-6D + 4A 
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Table 2: Energies in units of J obtained from approximate calculations on a quintet cluster 
of tetraliedra for AF exchange interactions. For details see the text. The degeneracy and 
spin of each state are indicated in parenthesis. In the row NS we indicate the number 
of total spin states and the number of local spin states present in each calculation. The 
total number of states with Mq = is 12870. The energies should be corrected by -12. 
The exact ground state energy has been calculated by |]59| . 



(5,D,Q,Mq) 


(0,0,0,0) 


(1,1,1,0) 


(1,1,1,0) 


(1,1,1,0) 


(1,1,1,0) 


exact 


Ssum 





1 


2 


3 


4 




NS = 


16/1296 


112/9232 


520/9232 


1192/9232 


1597/9232 


12870 


1 


-0.188(16;0) 


-0.235(16;0) 


-0.830(16;0) 


-0.830(16;0) 


-0.856(16;0) 


-0.82(16;0) 


2 




1.500(48;1) 


0.939(48;1) 


0.767(32;1) 


0.712(28,1) 


0.69(48;1) 


3 




1.812(32;1) 


1.690(32;1) 


0.938(16;1) 


0.750(1,1) 


1.23(32;1)) 


4 




2.797 (16;1) 


2.000(32;0) 


1.448(32;1) 


0.755(1,1) 


1.77(32;0) 


Est 




1.733 


1.768 


1.597 


1.568 


1.51 
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Table 3: The ground state energy per site and singlet triplet gap AEst, both in units of 
J. 



Eg 




-0.96 




-0.916 


1.688 


-1.00 




-1.144 




-0.975 




-1.0? 


1.10 


-1.12 




-1.10 


1.4 


-1.12? 


0.8 



method 



ref. 



J'/ J Expansion 
J' / J Expansion 
J'/ J Expansion 
Spin- Waves 
Mean Field 
J' / J Expansion 
J' / J Expansion 
Exact 16 site 
Super Tetrahedron: 16 sites 



5 

m 
m 
m 
m 
m 
m 
m 



Table 4: The b coefficients for s = 1/2 for the tetrahedron as basis cluster. 



s 





1 


1 


1 


1 


1 


1 





1 


1 


1 


2 


2 


2 


2 


2 


M 





-1 





1 


-1 





1 





-1 





1 


-2 


-1 





1 


2 


L 














1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


R 





1 


1 


1 











1 


1 


1 


1 


1 


1 


1 


1 


1 


m = -1/2 


1 











1/2 





-1/2 





1/4 





-1/4 


1/2 


1/4 





-1/4 


-1/2 


m = 1/2 


1 











-1/2 





1/2 





-1/4 





1/4 


-1/2 


-1/4 





1/4 


1/2 
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Table 5: Critical coordination number Zcrit and critical temperature in the pair-CVA 
Tc,pair/Tc^MF the tctrahedron-CVA Tc^h/Tc^mf for ferromagnetic nn Heisenberg exchange 
J and the ratio of the Ising and Heisenberg critical temperature T^j/Tc^h on the pyrochlore 
lattice - the Curie-Weiss mean field value is: fc^^c = 2JZs{s + l)/3 for Z = 6. The critical 
temperature in the Heisenberg limit (zero-field zero) {kBTc^n/ J) and in the Ising limit 
(ksTcj/J) for Dg/J = 5000. The ratio Tcj/Tc^h is in the mean field approximation this 
ratio is (1 + l/S')/3. The high temperature expansion is derived from the formula of 
Rushbrooke and Wood |Q . In the last two rows we give the characteristic critical values 
for the tetrahedral zero-field D^. 



s = 


1/2 1 3/2 2 5/2 3 7/2 


Pair-CVA with pairs, ferromagnetic Heisenberg spins: 


7 

Tc^pair / Tc^MF 


4 3 2.66 2.4 2.3 2.26 2.25 
0.607 0.737 0.773 0.788 0.796 0.801 0.804 


Tetrahedron- 


CVA with ferromagnetic Heisenberg (H) and Ising (I) spins: 


ksTcH/ J 
ksTci/ J 

Tc,H /Tc,MF 


-4.664 -9.636 -16.001 -23.980 -32.955 -43.566 
-8.808 -19.818 -35.234 -55.054 -79.277 -107.906 
0.583 0.642 0.667 0.679 0.687 0.692 
NA 1.888 2.056 2.208 2.295 2.406 2.477 


High temperature expansion Heisenberg model Z = 6 ferromagnetic: 


ksTj J 


0.629 0.684 0.699 0.705 0.709 0.711 0.712 
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Table 6: In this table we give the characteristic critical values for the tetrahedral zero-field 
Dr^max and the maximum critical temperature T^ax which characterizes the re-entrant 



re} 


2;iion for s = 1/2 to 7/2. 


See figure IC. 














Tmax/^(max) 


NA -2.647 


-5.128 


-8.143 


-12.400 


-15.662 


-20.137 




DT,niax/ J X 100 (max) 


NA -3.695 


-6.312 


-7.985 


-9.378 


-9.982 


-10.465 
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Figure 1: The pyrochlo (=B sublattice of the spinel structure) lattice. We have indicated 
the Corner-Sharing-Tetrahedra in blue, and some of the hexagon loops. The cubes (black 
thin lines) are drawn to guide the eye 

Figure 2: The tetrahedron with site and bond indices as used in the main text. The local 
trigonal (111) quantization axes and a tetragonal axis {C2 or S4) are also indicated. 

Figure 3: In this figure we illustrate the energy levels of spin states on a tetrahedron 
for s = 1/2. The exchange interaction is J = 1, the reduced tetragonal anisotropy is 
cIt = Dt/J, and the reduced exchange anisotropy is S = A/J. The reduced magnetic 
field is h = ^bH/J. The magnetic field is along a tetragonal axis. The {L,R) character 
of the states is indicated in parenthesis, the degeneracy is indicated in brackets. For the 
energies see also table: j^. 

Figure 4: The ground state phase diagram for ferromagnetic (F) and anti ferromagnetic 
(AF) nn exchange interactions for s = 1/2, in the (ItS plane. The numbers in parenthesis 
indicate the L, R character of the ground state. The equations separating the various 
ground states are given in the main text. 

Figure 5: The reduced critical temperature fc^Tc/ J as a function of Ct = 2{Zt — 1)/Zt for 
ferromagnetic nn exchange for s = 1/2; Zt is the coordination number. The case of 2, 3 
and 4 CST is Ct = 1, 4/3 and 1.5, respectively, are indicated as squares. 

Figure 6: Detail of the ferromagnetic phase diagram for s = 1/2 in the (ItS plane. We 
have indicate in black the region with re-entrant behaviour, which separates the singlet 
(S) ground state region from the ferromagnetic (F) ground state region. 

Figure 7: The reduced critical temperature ksTc/ J as a function of the reduced tetra- 
hedral anisotropy = Dt/J for various values of the reduced exchange anisotropy 
6 = A/J = 0., —1, —1.5, —2.0, —2.5, —3.0 for s = 1/2. Paramagnetic (P), Ferromagnetic 
(F) and singlet-like or collective paramagnetic (S) behaviour of the spin system is indi- 
cated. See also figure |^ and table ||. In this graph we have also indicated the curves, for 
the same set of 6 values, below which the entropy becomes negative: Tg. 

Figure 8: The number of states accessible as a function of the temperature for ferromag- 
netic exchange (F) and anti-ferromagnetic exchange (AF). 

Figure 9: The temperature T5 for which the entropy vanishes for non-magnetic spin 
system for ferromagnetic exchange J = —1 as a function of the tetragonal anisotropy dx, 
for magnetic fields h = 0, 1 and 2, along the line A-C in the phase diagram ^. 
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Figure 10: The reduced critical temperature ksTc/ J as a function of the reduced tetrago- 
nal anisotropy dx = Dt / J for s = 3/2 and reduced tetrahedral anisotropy. Paramagnetic 
(P), Ferromagnetic (F) and singlet-like or collective paramagnetic (S) behaviour of the 
spin system is indicated. The maximum temperature and anisotropy values given in 
table 1^. 
are indicated. 

Figure 11: The magnetization as a function of the magnetic field h for ferromagnetic 
exchange J = — 1 using equation p.46| , for three temperatures T = 0.1, 1 and 5, for various 
values of the tetrahedral anisotropy Irf^l = 0.05, 0.1, 0.25, 0.5, 1.0, 2.0 for 5 = —3 (line A-C 
in fig. ^ ) and 6 = 1 (line B-D in fig. g). A: 6 = —3, dx positive; B: 6 = 1, d^- negative; 
C: (5 = —3, dx'- negative; D:5 = 1, dx'- positive. For details see the text. 

Figure 12: The magnetization M. in the re-entrant region for ferromagnetic nn exchange 
J = — 1 for s = 1/2 in small field h = 0.001 for 6 = —3 for various values of the tetragonal 
anisotropy from dr = —0.85 to —1.075 with steps of —0.025. Units used: \J\. 

Figure 13: The magnetization as a function of temperature for ferromagnetic exchange 



J = — 1 using equation p.46|, for the magnetic fields h = H/\J\ = 0.01,1,2,3,5,10, for 
various values of the tetrahedral anisotropy dr = 0.05, 0.1, 0.25, 0.5, 1.0, 2.0 for 5 = 1 (line 
B-D in fig. ^). For details see the text. Units used: \J\. 

Figure 14: The specific heat as a function of temperature for small magnetic field in the 
re-entrant regiion for S = —3 and dr = 0.925. 
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